Entanglement from Charge Statistics: Exact Relations for Many-Body Systems 
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We present exact formulas for the entanglement and Renyi entropies generated at a quantum 
point contact (QPC) in terms of the statistics of charge fluctuations, which we illustrate with 
examples from both equilibrium and non-equilibrium transport. The formulas are also applicable to 
groundstate entanglement in systems described by non-interacting fermions in any dimension, which 
in one dimension includes the critical spin- 1/2 XX and Ising models where conformal field theory 
predictions for the entanglement and Renyi entropies are reproduced from the full counting statistics. 
These results may play a crucial role in the experimental detection of many-body entanglement in 
mesoscopic structures and cold atoms in optical lattices. 
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Introduction. — Entanglement entropy is playing an in- 
creasingly important role in describing quantum corre- 
lations in many-body systems Q. For a bipartite sys- 
tem, Fig. [l]a, the entanglement entropy of subsystem 
A is defined as Sa = — ^a{pa log pa}, where the re- 
duced density matrix pa = Ttb{/5} is obtained from the 
full density matrix p by tracing out the degrees of free- 
dom in the remainder B. For a pure state p = 
Sa = Sb = S. Entanglement entropy is currently being 
studied theoretically in a wide range of systems includ- 
ing quantum critical systems in one [Q-Q and higher |^| 
dimensions, topologically-ordered states g|, and evolu- 
tion after a quench [Q. Experimental progress, however, 
has so far been hindered by the difficulty of measuring 
the density matrix of a quantum many-body system and 
the fact that the definition of entanglement entropy itself 
does not refer to any directly measurable observables. 

There has thus been a growing interest in relating en- 
tanglement entropy to experimentally accessible quanti- 
ties, in particular fluctuations of charge and magnetiza- 
tion 0. An important step towards this goal was taken 
by Klich and Levitov who suggested an intimate con- 
nection between entanglement entropy and current fluc- 
tuations for non-interacting fermions M. They studied 
the entanglement entropy between two electronic leads 
connected via a quantum point contact (QPC), Fig. [l]b, 
and found that it can be expressed as a series in the 
cumulants of the current fluctuations. Remarkably, for 
purely gaussian fluctuations with variance Ci their result 
S = (7r 2 /3)C2 = (1/3) log(T/r) reproduces the confor- 
mal field theory (CFT) prediction S = (c/3) log(£/£) 
with the spatial extent C replaced by the temporal win- 
dow T during which the QPC is open, the spatial cut-off 
£ replaced by the short-time cut-off r, and the central 
charge equal to unity c = 1 |l0[ fflj . 

However, the series developed in Ref. ^ does not con- 
verge in general for non-gaussian fluctuations |l2j . This 
constitutes a major obstacle to further systematic study 
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FIG. 1: (color online), (a) The entanglement entropy charac- 
terizes non-local quantum correlations between subsystem A 
and the remainder B, such as (b) two electrodes upon opening 
of a quantum point contact (QPC) and (c) strongly-repulsive 
bosons hopping in a one-dimensional optical lattice. 



of entanglement and its relation to fluctuations in many- 
body systems. To illustrate the acuteness of this prob- 
lem, consider the simplest situation involving a QPC 
where only a single charge is transmitted with probability 
1/2: The series in Ref. ^expresses the entanglement en- 
tropy as an infinite sum over charge cumulants C n , where 
the coefficients for high-order cumulants asymptotically 
go to a constant (~ 2C n for large orders n). A quick 
calculation shows that for this simple process the C n di- 
verge as \C n \ ~ 2{n—l)\/ / K n , and indeed such factorial di- 
vergences are quite typical for non-gaussian distributions 
p3| . Experimental determination of the entanglement 
entropy from the cumulants using the results of Ref. ^ 
therefore requires further resummation. 

In this Letter we present exact and convergent series 
for the entanglement and Renyi entropies. The formu- 
las are applicable to a large class of many-body systems, 
which we illustrate with examples of entanglement gener- 
ation in a QPC, Fig. and groundstate entanglement of 
systems described by non-interacting fermions in one and 
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two dimensions. In one dimension the latter includes the 
equivalent, and experimentally relevant, case of strongly- 
repulsive bosons in an optical lattice, Fig. [l]c, as well as 
the spin-1/2 XX and Ising models. Our formulas pave the 
way to an improved understanding of fluctuations and 
entanglement in quantum many-body systems and may 
facilitate experimental investigations of entanglement en- 
tropy in mesoscopic structures where the measurement of 
high-order current correlators is becoming possible | Q |. 

Main result. — We begin by considering the QPC de- 
picted in Fig. Qb, but we will see below that the result 
applies to a much larger class of problems. Charge fluc- 
tuations in one of the leads of the QPC (the right lead, 
for concreteness) are characterized by the probability P n 
of transmitting n charges from the left to the right lead. 
The main result of this work relates the entanglement 
entropy to the charge statistics as 



iS = lim 



K+l 



t (K)C n 



(1) 



P n and x(A) = J2 n Pn& 1 is the generating function. 
The cutoff-dependent coefficients a n (K) are 



a n {K) 



'<)Sr K Si(k,n-1) 
L l~,k=n-\ k\k 



for n even, 
for n odd, 



(2) 



where Si(n, m) are the unsigned Stirling numbers of the 
first kind. A beautiful property of the series (0) is that 
only even-order cumulants contribute, reflecting the re- 
quirement that the entanglement entropy be symmetric 
between the left and right leads in the particular case 
of a pure state: Since charge conservation implies that 
the number of electrons n collected in the right lead is 
equal to —n charges collected in the left, only even-order 
cumulants are symmetric in the two leads. Moreover, 
a n (oo) = 2£(n) for even n where Ci n ) is the Riemann 
zeta function, and in particular we reproduce the CFT re- 
sult S = (7r 2 /3)C2 for purely gaussian fluctuations where 
all cumulants except the first and second are zero. In gen- 
eral, however, the number of cumulants included fixes the 
cutoff K . As K grows the series becomes an increasingly 
sharper estimate of S from below. This feature of the 
series is important from a practical point of view, since 
it ensures that including more terms always improves the 
estimate, e.g. Figs. || and ||. 

Derivation. — To derive Eq. ([!]) we begin with the ex- 
pression for the entanglement entropy 

S = -Tr {MlogM+ (1 - M)log(l - M)}, (3) 

where M = PrUjjPr is the correlation matrix projected 
onto the modes in the right lead by Pr, and we assume 
that M is either finite or can be regularized such that 
the operations below are well-defined jlBI. Here n is 



the Fermi-Dirac distribution before the evolution, and 
njj = UnW represents the state of the total system af- 
ter evolution U of the single particle modes. We first 
expand the logarithms in Eq. (||) around M = 0,1 to 
obtain S = Y] T A n /n with 



An 



Tr {M(l - M) n + M™(1 - M)}. 



(4) 



To relate the coefficients A„ to measurable quantities, we 
next use the Levitov-Lesovik determinant formula 
for the generating function of the charge transport statis- 
tics written as 



X(A) =dct[(l-M + Afe lA )e 



-i\UnP R U 



(5) 



For such a generalized binomial distribution it is 
useful to consider the factorial cumulants F n = 
9™ log x(— i log A)|x=i, since for n > 1 they are related 
in a simple way to Tr {M n } as F n = (-l)™" 1 ^ - 
l)!(Tr {M"}-Tr {UnP R W}). Solving for Tr {M n } and 
substituting into Eq. (^|) then gives 



where C„ = (— id\) n logx(A)|A=o are the cumulants of S = 



n=X 
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To write the sum in terms of ordinary cumulants, a cutoff 



K is introduced so that S = lim/f_ >00 X^=i A n /n. Using 
the relation F n = J2k=o(~^-) n ~ k Si( n > between fac- 
torial and ordinary cumulants, we finally arrive at Eqs. 
(|l|) and (0) after some algebra. 

Convergence. — The convergence of the series (|l|) for 
any M is most conveniently shown with the equivalent 
form ([}]). The counting statistics for non- interacting 
electron transport through a two-terminal conductor is 
always generalized binomial |L7|, such that the gener- 
ating function can be factored into a product of bino- 
mial events with individual probabilities < pi < 1. 
This is equivalent to evaluating Eq. (||) in the eigen- 
basis of M. Then S = J^i^iPi), where H 2 (x) = 
—xlogx — (1 — x) log(l — x) is the binary entropy func- 
tion. Since M is assumed to be finite and H2(x) has the 
convergent (absolutely, by the ratio test) series expansion 
££Li An/n with A n = Ei\Pi(l-Pi) n +Pi(l-Pi)l the 
series @ also converges to the correct value. Moreover, 
since each term A n /n is positive and decreasing for in- 
creasing n, Eq. ([3]) yields an increasingly sharper lower 
bound to the exact entanglement entropy. The same con- 
clusion holds for the series expressed in terms of ordinary 
cumulants, i.e., as the cutoff K is increased the sum (|l|) 
converges from below to the exact entanglement entropy. 

The rate of convergence generally depends on the dis- 
tribution of eigenvalues pi ; specifically, the expansion of 
the logarithms in Eq. (||) implies that the presence of 
eigenvalues of M near or 1 will slow the rate of con- 
vergence. As shown below, typically the first few terms 
dominate the series, but full convergence may require the 
inclusion of many terms. 
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FIG. 2: (color online). Exact entanglement entropy (dashed) 
and approximation by cumulants (solid) for a QPC with 
transmission D for (a) zero bias voltage and (b) bias volt- 
age V, scaled to the maximum value at D = 1 and D — 0.5, 
respectively. The cutoff K for the number of cumulants in- 
cluded, from bottom to top, is indicated in each plot. 



entropies. — Before turning to applications of 
we briefly show that the same approach can 
used to compute the Renyi entropies S„ = 
1 log[Tr {/S^}] of order n > 1, which can be 
as S n = (1 - n)- lr Tr {logfM" + (1 - M) n }}. 
similar derivation to Eq. (|l|) we obtain S n = 



^nR 



Pk(ji,R) 
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/ \mJ \nr — s/ (s— 1)\ ' 

for k odd. 

(7) 

Once again only even cumulants contribute to the sum. 
Moreover, for the gaussian case where all but the first and 
second cumulants vanish it can be shown that /^(n, oo) = 
(7r 2 /6)(l + 1/n), which correctly reproduces the CFT re- 
sult for the Renyi entropies Q . 

While we have expressed the entanglement and Renyi 
entropies in terms of the cumulants of the current with 
experimental applications in mind, it is also useful to ex- 
press the entropies in terms of the cumulant generating 
function itself. This can be done by relating both the 
Renyi entropies and the generating function to the spec- 
tral density fi(z) = Tr 8(z ~ M) of M to obtain 



du 



tanh u — tanh(nu) 



1 - n 



Im[logx(7r-2m)], 
(8) 

where the limit n — > 1 gives the entanglement entropy as 
1 



S = S 1 



du 



cosh 



-Im[log X (7T-2m)]. (9) 



By numerically integrating these expressions we can 
benchmark Eq. ([!]) against the exact results when only a 
finite number of cumulants are included, Fig. |^. 

Entanglement production in a QPC. — As our first ap- 
plication of Eq. (Q) we consider the generation of entan- 
glement due to current fluctuations in a QPC. For the 



equilibrium setup at zero bias and transmission D, the 
generating function is given by |19 



-AJG/(2tt 2 ) 



D sin 
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with G = \og{[h/3 / (ttt)] sinh[7ri/(/i/3)]}, where t is the 
time during which the QPC is open, r a short-time cut- 
off, and (3 = l/(fcsT) the inverse temperature. Here h 
and kj3 are the Planck and Boltzmann constants, respec- 
tively. For imperfect transmission D < 1, the generat- 
ing function ( |Toj ) possesses cumulants of all (even) or- 
ders, while for D = 1 the statistics is gaussian. At short 
times r«i< h/3 the growth of entanglement entropy 
is logarithmic, while at long times t h/3 the growth 
is linear. Thus at any non-zero temperature the entan- 
glement production rate eventually becomes constant. 
If a DC voltage bias V is applied, the non-equilibrium 
generating function at zero temperature (the quantum 
shot noise regime) is found to be the binomial distribu- 
tion logx(A) = (eVt/h)\og(l - D + De lX ), for which 
S = - (eVt/h)[D log D + (1 - £>)log(l - D)] @. This 
is also the electron-hole entanglement predicted for a bi- 
ased tunnel junction ppf , and indeed in both cases the 
Renyi entropies for small D are S n ~ [n/ (n — l)]N , where 
N = (eVt/h)D is the number of transferred particles [^lj . 

The dependence of the entanglement entropy on the 
transmission D in the two cases is shown in Fig. ||. We 
see that the first few cumulants give a good approxima- 
tion to the exact result obtained from Eq. (||), especially 
for the case of a bias voltage. This indicates that the 
detection of entanglement generation in a QPC may be 
within experimental reach using currently available noise 
measurement techniques p4| . 

Groundstate entanglement in the XY model. — We next 
consider a very different setup from the QPC, namely 
the bipartite entanglement entropy of the groundstatcs 
of systems described by free fermions, such as the spin- 
1/2 XY model in one dimension with Hamiltonian 



H = J^[(l+7)5f^ 1 + (l- 7 )Sf% 1 ] 



(ii) 



This Hamiltonian is equivalent to a problem of free spin- 
less fermions through the Jordan- Wigncr transformation 
p2[ and for 7 = also describes bosons in a one- 
dimensional lattice with infinite on-site repulsion, Fig. 
|l|c. Here we are concerned with relating the entangle- 
ment entropy to the fluctuations as a function of subsys- 
tem size I. The central quantity is M = (1 + v / G T G)/2, 
where Gij = ((aj — ai)(a] + hj)) with the indices i,j 
restricted to the first £ sites of the chain and di is the 
Jordan- Wigner fermion annihilation operator on site i. 
The entanglement entropy is then given by Eq. (||) M, 
while the corresponding generating function is x(A,^0 = 
det(l — M + Me lX ) up to an irrelevant phase. For the 
XX model at 7 ~ the generating function describes 
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FIG. 3: (color online). The exact entanglement entropy and 
approximation by cumulants for (a) the critical spin-1/2 XX 
chain (upper set of curves) and critical Ising chain in trans- 
verse field (lower set of curves) with periodic boundary condi- 
tions, L = 100, and (b) free fermions on L x L square lattices 
with periodic boundary conditions, with the partition taken 
to be half the lattice. Solid curves are the exact entangle- 
ment entropy, while the cutoff number K for both figures is 
2 (circles), 4 (squares), 8 (crosses), and 40 (plusses). 

fluctuations of S z , i.e., x(A, i) = (exp(iA J2i=i $f )) = 
(exp(iA ^ i=1 (ajaj — 1/2))), while for 7 7^ the fluc- 
tuations must be interpreted as those of quasi-particlcs, 
since total S z is not a conserved quantity. The results 
for both the critical XX model at zero magnetic field 
(7 = 0,6 = 0) and critical Ising model (7 = 1,6 = J) 
with periodic boundary conditions are shown in Fig. |^a. 
Remarkably, the cumulants correctly reproduce the two 
central charges c = 1 and c = 1/2 for the XX and Ising 
models, respectively We also note that the random 
singlet phase || can be studied with this method by ran- 
domly varying J and 6 across different sites. 

More generally, the same formalism is applicable to 
any quadratic fermionic Hamiltonian, regardless of di- 
mension, and as such can also be used to determine the 
entanglement entropy from the fluctuations in higher di- 
mensions. On the L x L square lattice in two dimensions 
with nearest-neighbor hopping, for example, we find that 
the approximation by cumulants follows the predicted 
S ~ L log L scaling |2^| as shown in Fig. ||b with the 
subsystem taken to be half the lattice. 

Conclusions. — We have derived exact formulas for the 
entanglement and Renyi entropies in terms of the statis- 



tics of charge fluctuations. The expressions are applicable 
to a wide range of experimentally relevant many-body 
systems, for example entanglement produced in QPCs 
and groundstatc entanglement of systems described by 
non-interacting fermions (or quasi-particles) in any di- 
mension. 
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